tmax = 2000; l = 32;
alpha = .08; beta = .0; gamma = .05; mu_v = .005; mu_is = .005;
Q_v = (1 - mu_v)^l;
Q_is = (1 - mu_is)^l;

